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We investigate the yield ol A(1520) resonance in heavy ion collisions within the framework of a 
kinetic master equation without the assumption of chemical equilibrium. We show that reactions 
such as A(1520)+7r ^ E* can favor E* production, thereby decreasing the A(1520) yield. Within the 
same approach we thus find a yield enhancement for E(1385) and a yield suppression for A(1520). 
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I. INTRODUCTION 

Hadron resonances are observed in a surprisingly large 
yield when a quark-gluonplasma (QGP) fireball breaks 
up into hadrons [H, 0,11, 0, H, @] ■ This is unexpected, since 
the invariant mass signature formed from decay products 
could be erased by rescattering of the strongly interact- 
ing decay products 0|- Thus a direct detection of res- 
onances implies an exceedingly short period of hadron 
scattering, and/or a final state repopulated by hadronic 
interactions Q. As a result, the final resonance yield 
can be considerably different from statistical hadron gas 
(SHG) benchmark expectation. It has already been re- 
ported that the short lived (compared to characteristic 
hadron phase evolution times) resonances are in general 
enhanced Q compared to SHG benchmark yield. 

The new result, we obtain here, is that the long lived 
resonances, such as A(1520), can be considerably sup- 
pressed in their yield. This effect is amplified for the 
case when the initial hadron fugacities, and thus particle 
yields, are above chemical equilibrium. This situation is 
expected for a hadronizing QGP phase. The low A(1520) 
yield has been reported both in RHIC and SPS experi- 
ments [H, 

In a global QGP breakup (hadron chemical freeze-out) 
particles and resonances are formed. Many resonances 
have a relatively large decay rates. This implies a large 
scattering formation rate. Thus resonances are excep- 
tionally strongly interacting particles and continue to 
evolve in what we call kinetic phase, even after all other 
particles freeze-out. This continued reaction phase is spe- 
cific to the resonances and can last considerably beyond 
the last non-resonant (elastic) scattering. 

The resonance suppression, or enhancement, mecha- 
nism works as follows. In thermal hadronic gas the reac- 
tion, 

1 + 2^3, (1) 

can occur in both directions: the resonance decay 3 
1-1-2, and the back-reaction (regeneration) resonance for- 
mation 1 + 2 3. When the reaction goes with the 
same rate in both directions, we have chemical detailed 



balance, e.g. particles yields do not change in this pe- 
riod of temporal evolution of the system. This does not 
necessarily mean that we have a chemical equilibrium. 
Instead it may be a transient condition for which none 
of the three particles is equilibrated chemically - we will 
show when this can happen. 

In the study of resonance decay and regeneration we 
are using the momentum integrated population master 
equations. We assume a fireball expansion model gov- 
erned by hydrodynamic inspired flow with conserved en- 
tropy content. In our considerations we presume that 
the yield of pions it is so large that we can assume it 
not to be materially affected by any of the reactions we 
consider. Thus we fix pion yield in terms of an ambi- 
ent fugacity and temperature value, and in essence the 
total (per unit rapidity at RHIC) yield is fixed since we 
conserve entropy. 

An important assumption implied below is that the 
rapidly expanding hadron system maintains for the rel- 
evant particles a fully thermal (Boltzmann) momentum 
distribution. To describe the evolution of hadron abun- 
dances in the kinetic phase we track in time the yields 
of single strange hadrons after their initial formation. 
This is implemented in terms of time dependence of the 
chemical fugacities T(t), and the time dependence of the 
hadronization temperature T{t). 

We look in detail at three potential evolution scenar- 
ios: 

a) a high temperature breakup at Tq ~ 180 MeV where 
the entropy content of the equilibrated QGP and HG- 
phase are similar; 

b) the To ~ 160 MeV case where chemical non- 
equilibrium among produced hadrons is already required; 
and 

c) at To ~ 140 MeV which is favored by descriptions 
of stable hadron production, and in which case a strong 
chemical non-equilibrium situation arises. 

For the late stage of the expansion, at relatively low 
density the assumption of thermal momentum distribu- 
tion may not be anymore fully satisfied. In particular 
pions of high momentum could be escaping from the fire- 
ball. For this reason we will consider here a second sce- 
nario, which we call "dead channel" . In this scenario we 



2 




FIG. 1: (color on line) Reactions scheme for A(1520) and E(1385) population evolutions. 



assume that the reaction ([T]) goes mainly in the direc- 
tion of resonance 3 decay and the resonance formation is 
switched off for 

ma - (mi + m2) > 300 MeV. (2) 

Without a complete kinetic model including equilibra- 
tion and particle emission we do not know the exact en- 
ergy in condition ^ and timescale (during expansion) 
for which Boltzmann distribution is violated and dead 
channels appear. It is possible that reality lies between 
the two cases (kinetic Boltzmann distribution and dead- 
channels) considered here which, in our opinion, are the 
two most extreme limits. 

For A(1520)/A'^ ratio calculations we employ and de- 
velop further the approach used for I](1385)/A'' in Q. 
However, in chapter III Al we investigate many further 
reactions in which resonance A(1520) participates. Thus 
we are obliged to develop a completely numerical evolu- 
tion, for which the analytical study of S(1385)/A° pro- 
vides a benchmark check of our approach. In addition 
to new and numerous reaction channels we also intro- 
duce deformation of the reaction rates due to stimu- 
lated Bose enhancement of the reactions. This formal- 
ism is presented in chapter III Bl We discuss the tem- 
poral evolution of HG particle fugacities T(i) in IIII Al 
In chapter IIIIBI we present results for the evolution of 



particle S(1385), A(1520) multiplicities during kinetic 
phase. In chapter IIII CI we obtain the observable 'ob' 
ratios A(1520)ob/Atot and I](1385)ob/Atot- We discuss 
our results in section HVl 

II. KINETIC EQUATIONS 

A. Reactions scheme for A(1520) and E(1385) 

In figure [1] we show the scheme of reactions which all 
have a noticeable effect on A(1520) yield after the chem- 
ical freeze-out kinetic phase. The format of this presen- 
tation is inspired by nuclear reactions schemes. On the 
vertical axis the energy scale is shown in MeV. There 
are three classes of particle states, which we denote from 
left to right as "iV" (S=0 baryon), "S" (S* = -1, / = 1 
hyperon) and "A" (S* = — 1, / = hyperon). Near each 
particle bar we state (on-line in blue) its mass, and/or an- 
gular momentum and/or total width in MeV. The states 
A(1520) and S(1385) are shown along with the location 
in energy of A(1520) -I- tt and E(1385) + tt respectively, 
both entries are connected by the curly bracket, and are 
highlighted (on-line in red) . The inclusion of the 7r-mass 
is helping to see the kinetic threshold energy of a reaction. 
The lines connecting the iV, E, A columns are indicating 
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the reactions we consider in the numerical computations. 
All reactions shown in figure [T] can go in both directions, 
as shown by the double arrows placed next to the numer- 
ical value of the partial decay width Ti in MeV. 

A(1520) decays with a total decay width of about 15.6 
MeV, with two main channels: 



E- 
7V- 



TT ^ A(1520), 
K ^ A(1520), 



r f 
r 



6.5 MeV; 
^ 7 MeV. 



(3) 



However, A(1520) reacts with several heavier S*- 
resonances, (E* = E(1670), E(1750), £(1775), E(1940), 
S(2030)): 



A(1520) + 7r^E*, 



(4) 



and these reactions have a larger reaction strength 
shown in figure [TJ A(1520) nearly behaves like a 'stable' 
hadronic particle since: 

a) it is dominantly coupled to heavier resonances; 

b) its natural lifespan is larger than the hadronic 
reaction rate. 

Hereto we note that (several) S* involved in Eq. ^ 
participate in further reactions: 



A(1115) + 7r ^ E*; 


(5) 


S(1190) + 7r ^ E*; 


(6) 


N + K ^T.*; 


(7) 


E(1385) + 7r ^ E*; 


(8) 


A + /s:^E(1940,2030); 


(9) 


N + K{892) ^ E(1940); 


(10) 


E + r? ^ E(1750). 


(11) 



All reactions shown above can excite E* resonances. 
Since the mass of A(1520) is near to the E* mass, the 
yield of A(1520) is effectively depleted by the reaction 
chain 



A(1520) 



7V + K, etc. 



(12) 



The balancing two step back-reaction can also occur, es- 
pecially once A(1520) has been depopulated. Thus a dy- 
namical reduced detailed balance yield of A(1520) would 
result if the system were at fixed volume rather than ex- 
panding. 

The multiplicity of E(1385) is mostly determined by 
its dominant decay and production in the reaction 



A(1115) + 7r ^ E(1385), 
and to a lesser extent by the reaction 

E(1190) -Htt ^ E(1385). 



(13) 



(14) 



The resonance E(1385) participates further in reactions 
with heavier E*; see reaction (|8]), but strength of these 
interactions is smaller than for similar reactions with 
A(1520) and smaller than the decay width of E(1385). 



Thus we find that the influence of these reactions on 
E(1385) yield is small. Another reason for a reduced ef- 
fective depletion rate of E(1385) is that a lesser fraction 
of this resonance is needed to excite E*. Thus in such 
a reaction the depopulation effect decreases because of a 
larger mass difference between S(1385) and E* in com- 
parison with A(1520) and E*. 

The reactions scheme for A(1520) reactions with dead 
channels is shown in figure [2l The difference between 
figure[T]and figure[2]is that some of the reaction lines have 
single-directional arrows, as is stipulated by the condition 



B. Resonances densities, time evolution equations 

The evolution in time of the resonance yield is de- 
scribed by a master equation, where the process of reso- 
nance formation in scattering is balanced by the natural 
resonance decay: 



1 dN3 

v~dr 



^ dVdt 



■1+2 



dVdt 



(15) 



where subscripts i, j denote different reactions channels 
when available. We further allow different subscripts 
i, j for the case where there are dead channels. Thus 
dWl^2^3/d-^'^t and dWl^i^2/'^^d,t are invariant rates 
(per unit volume and time) for particle 3 production and 
decay respectively. In case all reactions occur in both di- 
rections the total number of fusion channels is the same 
as the total number of decay channels. 

Allowing for Fermi-blocking and Bose enhancement in 
the final state, where by designation particles 1 and 3 are 
fermions (heavy baryons) and particle 2 is a boson (often 
light pion) we have for the two rates: 



dWi_ 



•1+2 



dVdt 



-h 



d^Pi 



d^P2 



2Ei{2'k)' 



(l-/i) 



2E2 (27r) 



:-5]|(p3|M^'| 

53 ^ 



PlP2) 



(16) 



J 2Ei{2nY^^J ' 



52rf^P2 



r/2 



dVdt J 2Ei{2TTf-'"J 2E2{2iTf- 

(1 - /3) (2^)^ 5' ip, +P2- P3) 



2E3 (2^)^ 

— J2\{p,P2\M' 



9i92 



P3) 



(17) 



spin 



where gi,i — 1,2,3 is particles degeneracy. The Bose 
distribution function for particle 2 is 

= T-vi/^-i' ^^^^ 
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FIG. 2: (color on line) Reactions scheme for A(1520) and S(1385) interactions in the "dead channel" model. 



and Fermi for particles 1, 3 are: 



1 



J = 1,3. 



(19) 



Here T,; is particles fugacity, and u ■ pi = Ei, for 
= (1,0) in the rest frame of the heat bath where 
d'^p5o{pj — rnj) (Ppi/Ei for each particle. Hence, 
Eg. p^ and Eq. (fT7)l are Lorentz invariant, and thus as 
presented these rates can be evaluated in any convenient 
frame of reference. Normally, this is the frame co-moving 
with the thermal volume element. 

For the heavy baryon (resonances), particles 2,3, we 
can work using the expansion of the relativistic distribu- 
tion, the first term is the Boltzmann limit: 



(20) 



where Xi = rrii/T, K2{x) is the Bessel function (not to be 
mixed up with particle 2) . However, we use the complete 
Bose distribution to describe pious. 

We introduce in medium lifespan of particle 3: 



1 

T3 



123 



V-^dNs/dTa 



(21) 



and, similarly, channel lifespan T3, omitting the sum 
Here the rate i?i23 is: 



R 



123 



d^pid^P2d^P3 



8EiE2E3{2t:) 



-1 
3 



X 5^ {pi +P2-P3) e"-P^/^ I ^P^P^ \M'\p3)f- (22) 

spin 

R is independent of the fugacity in the Boltzmann- limit. 
In next section III CI we will see in what way the in- 
medium decay rate varies from the free space decay rate. 
This is due to the effect of quantum enhancement, and 
the fact that a particle emerged in a thermal bath with a 
finite temperature. A particle is not decaying in its rest 
frame. 

The production and decay rates are connected to each 
other by the detailed balance relation [ll|, [l2| : 



-i^-idWi+2- 



dVdt 



= T 



-idW3 



dVdt 



= i?i 



23- 



(23) 



Using detailed balance Eg. (1^51) we obtain for fugacity T3 



m. pap 1 

the evolution equation [ll|, |l2| : 

dT ^ Tl \tt ts ^rij 



(24) 
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where we have also introduced characteristic time con- 
stants of temperature T and entropy S evolution 



1 

TT 

1 

TS 



dT 

d\n{VT^) 



dT 



T. 



(25) 
(26) 



The entropy term is negligible, ts S> T3, tt since we im- 
plement near conservation of entropy. We implement this 
in the way which would be exact for massless particles 
taking VT^ =Const.. Thus there is some entropy growth 
in HG evolution to consider, but it is not significant. In 
order to evaluate the magnitude of tt we use the re- 
lation between Bessel functions of order 1 and 2 (not 
to be mixed up with particles 1,2) d[z^K2{z)) /dz = 
—z'^Ki{z). We obtain 



TT 



(27) 



Tt > 0- We invoke a model of matter expansion of the 
type used e.g. in [lB|, where the longitudinal and trans- 
verse expansion is considered to be (nearly) independent. 
In this model we have: 



T 
T 



1 f2{vT/Rj 
3 



(28) 



where R± is the transverse radius, v is the velocity of ex- 
pansion in the transverse dimension. All flow parameters 
(or temperature dependence on r) are the same as in 
For a static system with — > we see that Eq. ([M)) has 
transient stable population points whenever 



(29) 



Next we address the functional dependence on time 
of Ti,T2. In the equation for Ti we have terms which 
compensate what is lost/gained in T3 see Eq. (p4)) . Fur- 
ther we have to allow that particle '1' itself plays the 
role of particle 3 (for example this is clearly the case for 
A(1520)). That allows a chain of populations relations 
as follows: 



{!' + 2' ^l) + 2^ 3, 



Then we obtain: 



dTi 



^ 1 dN§/dTl 1 



(30) 



„ 1 dN^/dT^ 



Ti 



[tt ts ^Tij ^ ' 'rl 



(31) 



of equations for baryons closes with the equation for Ti 



dTy 



Ti 



1 dN^/dT'l 
^ dNv/dTv 



1 

Tt 



1 

Ts 



tI' dNv/dTv 
(32) 



In the present setting T2=7r =Const. by virtue of entropy 
conservation (see discussion below) and the same applies 
to the case 2' = tt. However, if either particle 2 or 2' is a 
kaon, we need to follow the equation for T2.2'=/f which 
is analogous to equation for particle 1 or 1'. 

The evolution equations can be integrated once we de- 
termine the initial values of particle densities (fugacities) 
established at hadronization/chemical freeze-out. We 
determine these for RHIC head-on Au-Au collisions at 
Y^snn = 200 GeV. We introduce the initial hadron yields 
inspired by a picture of a rapid hadronization of QGP 
in which quarks combine into final state hadrons. For 
simplicity we assume here that the net baryon yield at 
central rapidity is negligible. Thus the baryon-chemical 
and strangeness potentials vanish. The initial yields 
of mesons (gg, sq) and baryons (ggg, qqs) are controlled 
aside of the ambient temperature T by the constituent 
light quark fugacity 7, and the strange quark fugacity 7^ . 

The strangeness pair-yield in QGP is maintained in 
transition to HG. This fixes the initial value of 7s. In 
fact, since we investigate here relative chemical equilib- 
rium reactions our results do not depend significantly 
on the exact initial value 7s and/or strangeness content. 
The entropy conservation at hadronization fixes 7^. For 
hadronization temperature T{t = Q) = Tq ~ 180 MeV, 
7q = 1. However, when Tq < 180 MeV, 7, > 1 in or- 
der to have entropy conserved at chemical freeze-out. At 
To — 140 MeV 7^ = 1.6 that is close to maximum pos- 
sible value of 7q, defined by Bose-Einstcin condensation 
condition p^ . 

For reactions, such as shown in Eq. ([1]), we have (lower 
index defines particle considered, where y = E, A is a 
hyperon) 



JL(i=y) - 



ills 



(2=7r) 



lo 



^ (2=K) 



Iqls 



(33) 



(34) 



where the particle 1 in reaction ([T]) is a baryon and par- 
ticle 2 is a meson. The particle 3 is always a strange 
baryon: 



-Y-o 

(3=1-) 



ills 



Note that for 7^ > 1 we have always initially 

7a > 1- 



(35) 



(36) 



The ratios of derivative of Ni seen in the first line are due 
to the definition of relaxation time Eq. (HI]). The system 



As a consequence initially the pair of particles 1,2 reacts 
into 3. 
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As already noted, we do not need to follow the evo- 
lution in time for the pion yield, which is fixed by con- 
servation of entropy per unit rapidity, as incorporated in 
Eq. (|28p. Thus it is (approximately) a constant of mo- 
tion. This can be seen recalling that the entropy per 
pion is nearly 4 within the domain of temperatures con- 
sidered. Thus the conservation of entropy implies that 
pion number is conserved. With VT^ ~ Const., this fur- 
ther implies that during the expansion 

= 7g = Const., 

which we keep at the initial value. 

C. In medium lifespan calculations 

In our calculations we take into account the influence 
of the medium on resonance lifespan and the effect of the 
motion of the decaying par ticle with respect to the ther- 
mal rest frame. In [ll|, [Tj] it was noted that the decay 
rate i?i23, Eq. ([22]), of particle 3 (density) in a thermally 
equilibrated system can be cast into a form which in- 
volves the free space decay rate: 



-0-123 — 



To Jo 



1>(P3), (37) 



where function ^{ps) for reaction (fT3| 2(1385) ^ A + tt 



IS 



^(Ps) 



;TAe''-|-e-°^)(e°i-T^e- 



In 



ai 



ElE^ 



02 



6 = 



V P3 



(38) 



Here p* — pi — p2 and ^ 



* 2 



are the mag- 



nitude of the momentum and, respectively, the energy, of 
particles 1 and 2 in the rest frame of the particle 3. From 
energy conservation: 



E- 
V 



m| ± (rn\ — m^) 



1,2 
* 2 



2m3 



- Eu2 



'1,2 



{ml 



Ami 



(39) 



Here fugacities for A and tt correspond to those for 
particles 1 and 2, respectively. For the temperatures 
of interest (hadronization of QGP and below) m\ and 
ms » T . With sufficient accuracy we can write 



1 



In- 



-L -7] 



-6 



(40) 



There are no significant medium effects upon decay rate 
of I](1385) and A resonances. However the pions have 
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FIG. 3: (color on line) The Boss enhancement factor 1 -|- 
f-n{E1) in E(1385) rest frame as a function of light quark 
fugacity 7, for the reaction E(1385) ^ Att at T = 140 MeV 
(blue, solid line), at 160 MeV (green, dash-dot line) and 180 
MeV (red, dashed line). The dots show the initial value of 
fugacities for the three possible hadronization cases. 



energy E^ = 250 MeV (Eq.®) in the S rest frame and 
the Bose enhancement effect is possible in the oversatu- 
rated hadronic gas after QGP hadronization. 

For the low temperatures considered here we can as- 
sume that E resonances almost do not move. Thus the 
enhancement effect in the thermal bath frame is close to 
the enhancement in the 2(1385) rest frame. The decay 
rate increases by Bose enhancement factor 1 + f-^ (here 
/jr = f^{E2, T)). In figure [3] we show Bose enhancement 
factor as a function of light quark fugacity 7^ for tem- 
perature To ^ 140 MeV (blue, soHd Hue), Tq = 160 MeV 
(green, dash-dot line). To = 180 MeV (red, dashed line). 
The large dots show Bose enhancement factor for our 
initial 7^ determined from entropy conservation in fast 
hadronization. The fugacity 7^ = 1.6 is close to maxi- 
mum expected value at To = 140 MeV. The maximum 
fugacities for each temperature correspond to Bose - Ein- 
stein singularity. The Bose enhancement effect is largest 
for maximum 7^ and it diminishes for small 7^. At fixed 
entropy the greatest enhancement is for smallest ambient 
temperature, see the dot on solid line in figure [31 

In figure [4] we show the corresponding decrease in the 
lifespan, the ratio t^/tq as a function of temperature T 
in the reaction E(1385) ^ Att. We consider temperature 
range from corresponding hadronization temperature un- 
til T = 70 MeV. We assumed, that Tjr is a constant. 
Fugacities of heavy resonances do not influence the re- 
sult. The lowest ts/tq ratio is for 7^ — 1.6 at To ~ 140 
MeV when we have maximum value of 7^ for given tem- 
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80 100 120 140 160 180 
T [MeV] 

FIG. 4: (color on line) The ratio of the in medium lifespan 
rs with the vacuum lifespan tq as a function of temperature 
T for the reaction E(1385) « Att. The dashed (red) line is 
for hadronization at To = 180 MeV, 7, = 1.0; the dot-dashed 
line (green) for hadronization at 160 MeV, 7q = 1.27; solid 
line (blue) is for hadronization at 140 MeV and jq — 1.6. 



perature. If we compare this value of t^/tq = 0.65 with 
inverse Bose enhancement factor l/{l+fT^{E2,T)) — 0.54 
for this T and jq (see figure [3]) we see that these values 
are near to each other (difference is about 20% ) as ex- 
pected for ms >> T. For smaller T, decay time goes 
to its vacuum value. 

The same calculations are applicable for heavier S*. 
When the difference of mass of the initial and final state 
resonance decreases, the Bose enhancement effect in- 
creases, since it involves small momenta. The largest 
effect is for reaction i;(1670) ^ A(1520) + tt. On the 
other hand, for the reactions which satisfy condition ^ 
the enhancement effect becomes so small that we do not 
need to include it in our calculations. 



solve this system of equations numerically, using classical 
fourth order Runge-Kutta method. 

Particle fugacities, except T^r, change rather rapidly. 
Figure [5] shows the computed T{t) as a function of tem- 
perature T(t). We present here the scenario in which all 
reactions evolve in both directions, for the initial condi- 
tion js — Jq- The time, corresponding to the tempera- 
ture shown at the bottom, is shown at the top of figure [5l 
in each frame. On the left we have hadronization at 140 
MeV, in the middle at 160 and to the right at 180 MeV. 
Each frame has the same scale size for temperature unit, 
not time. For Tj;. we show two possible evolution exam- 
ples, for Si75o (dash-dot dark fine) and S(1775) (dashed 
line). These resonances have significant influence on the 
A(1520) yield. The sohd lines are for Ta(i520) (upper, red 
line) and Ts(i385) (lower, hght blue hue). The dash-dot 
and dashed light lines are for T^j-iigQ) and T^o, respec- 
tively. The upper dotted line is for T^v and lower dotted 
line is for Tk- 

An important feature is that the Ts of massive hadron 
(resonances) increase very fast when T decreases. This is 
so since in absence of a rapid re-equilibration reactions, 
multiplicity of given resonance must be conserved. Then, 
according to Eq. pp]) oc 1/K2{mi/T), and thus for 
large rrii cx exp(mi/r). We would expect T; > Tj, 
when TTii > rrij, and T decreases. This behavior is jus t 
like we found for the case of large charm fugacity [10| . 
However, because of the decay and regeneration reac- 
tions there are some deviations from this expectation in 
figure O 

For To = 180 MeV in most cases T3 > T1T2 (i > 0). 
Massive resonances decay to lower mass particles. The 
result is defined by resonance mass, its decay width 
and decay products. For example Tx;(i775) is smaller 
than T5](i75o) and Ta(i52o), because of its large decay 
width. Therefore excitation of S(1775) by A slightly 
dominates over I](1775) decay to A(1520) even in this 
case, when for most resonances the decay is dominant. 
For smaller initial hadronization temperatures Ta(i52o) 
becomes smaller than and even smaller than 

T"e(i385) in some range of temperatures. This suppres- 
sion occurs because of I](1775), and others S* regener- 
ation. Because of large T^, Tj](i775) < ^^(1520)^,^ the 
E(1775) production by A(1520) is dominant in the full 
range of T considered here. 



III. NUMERICAL RESULTS 

A. Evolution of fugacities 

In order to evaluate the A(1520) and £(1385) mul- 
tiplicities we must integrate Eq. (fM)) . or Eq. (PT|) . or 
Eq. (|32p for each particle involved in figure[Tl and perform 
similar operations for reactions with dead channels in fig- 
ure [3 This system of equations includes equations for 
A(1520), E(1385), five equations for I]*s, equations for 
K(892) and A and equations for ground states A(1115), 
S(1190), N, K. All reactions in figures[T]are included. We 



B. Final A(1520) and E(1385) multiplicities 

In this section we consider the evolution of the multi- 
plicity of resonances A(1520), £(1385), 11(1775) during 
the kinetic phase. We use the Boltzmann yield limit, 
Eq. (|20|) . By the symbol X{T) we refer to a particular 
resonance, and Xq is the initial multiplicity for that res- 
onance. The dynamic yield of this resonance may be 
expressed as 

XjT) ^ Tx{mt)^K2{mxint)) 
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FIG. 5: (color on line) The fugacities T for selected particles are shown as a function of temperature T{t), for To = 140 MeV 
on the left, for To — 160 MeV in the middle and for Tb = 180 MeV, on the right. See text for further details. 




FIG. 6: The ratio E(1385)/S(1385)o on left and A(1520)/A(1520)o on right as a functions of temperature T{t) for different 
initial hadronization temperatures To — 140, 160 and 180 MeV (blue/bottom, black/middle and red/top lines, respectively). 
Solid lines are for calculations with dead channels, dashed lines are for calculations without dead channels. 
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FIG. 7: The rates for main channels of E(1775) (on the left) and E(1750) (on the right) decay and production as a functions 
of temperature T in the case when all reactions go in both directions and To — 140 MeV. Solid lines are for reaction E* « 
A(1520) + tt; dash-dot lines are for reaction E* ^ N + K; dashed lines are for reaction E(1775) <-> A" + vr on the left and 
E(1750) ^ E + on the right; blue and red lines are for decay and backward fusion reaction, respectively. 



Figure [S] shows this yield as a function of T(t) for X — 
E(1385) (left) and X = A(1520) (right). We consider 
three initial conditions, temperature Tq = 140, 160, 180 
MeV, with corresponding 7^ = 1.6, 1.27, 1.0, respectively. 
The solid lines correspond for the model with dead chan- 
nels and dashed one are for case when all reactions are 
symmetric in both directions. The thin dotted vertical 
line at T = 120 MeV marks the kinetic freeze-out temper- 
ature, assumed before in i9]. The main result is that the 
resulting relative yields for A(1520) and 2(1385) behave 
qualitatively different from each other. In particular, as 
the temperature decreases, for the case Tq = 140 MeV 
we observe a strong yield suppression for A(1520), and a 
strong enhancement for I](1385) (as compared to initial 
SHM yields). 

To better understand the mechanism of A(1520) sup- 
pression, we analyze in some detail the case of S(1775) 
and S(1750) decay and production rates dW/dVdt. We 
assume here that these reactions can go in both direc- 
tions. In figure [7] we show the reactions rates for the 
principal channels of decay and production as a functions 
of temperature T for E(1775) (left) and S(1750) (right), 
for the case of initial temperature Tq = 140 MeV which 
provides the largest A(1520) suppression. Solid lines are 
for the reaction S A(1520) -f tt, dash-dot lines are 
for reaction Yj N + K , dashed lines are for reaction 
S ^ A" + TT. Two set of lines are presented for the decay 



(on-line blue) and backward fusion reaction (on-line red), 
respectively. 

As temperature decreases, all rates dW/dtdV are in- 
creasing rapidly. This is mainly because fugacities T 
increase nearly exponentially when number of particles 
is conserved, see figure [S] We see that at the begin- 
ning of the kinetic phase all reactions go in the direction 
of S(1775) production, since S](1775) production rate is 
larger than its decay rate for all channels. Then at first 
S(1775) ^ + TT decay rate becomes dominant over 
S(1775) production rate in this channel, followed by the 
same for E(1775) ^ N + K channel. 

For the reaction i;(1775) ^ A(1520) -I- tt backward re- 
action is always dominant. As result, during the kinetic 
phase always more A(1520) resonances are excited into 
i;(1775) than they are produced by E(1775) decay. The 
reason for this is the decay of S(1775) to the other chan- 
nels, as long as Ts(i775) < TI"a(i520)T^. The hghter is 
the total mass of decay products, the earlier the decay 
reaction becomes dominant. This is due to the fact that 
the fugacity of T for heavier particles increases faster 
with expansion. Therefore, the decay rate becomes dom- 
inant earlier, when the difference between initial and final 
mass is larger. The net result is A(1520) suppression by 
i;(1775) excitation. 

In figure [5] we show the yield of 2(1775) normalized 
by its initial yield at hadronization: E(1775)/S(1775)o 
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as a function of T{t). Like in the other figures above, 
sohd lines are for the dead channels and dashed lines 
are for case when reactions go in both directions, solid 
(blue) lines are for Tq = 140 MeV, solid (black) lines 
for To = 160 MeV, and sohd (red) lines are for Tq = 180 
MeV. Each of the lines can be identified by their initial T- 
value. We see that when all reactions go in both direction 
the ratio E(1775)/I](1775)o increases at first similar to 
I](1385)/S(1385)o and A(1230)/A(1230)o ratios [|. 

Compared to these ratios, I](1775)/I](1775)o ratio 
reaches its maximum value earlier, and after the max- 
imum, the yield of 2(1775) decreases faster. The reason 
for this behavior is that the mass of E(1775) is larger. 
The phase space occupancy ^^.(177^)1 ^-nd therefore its 
decay rates, increase faster than the fugacity and de- 
cay rates for S(1385) and A(1230). Therefore decays 
S(1775) to some channels and its total decay rate be- 
come dominant earlier (see figure [7]). Although the total 
decay width of E(1775) is approximately the same as for 
A(123Q), the maximum value of this ratio is smaller. 

Said differently, the maximum yield of S(1775) does 
not have time to reach the value as high as that for 
A(1230). We thus learn that the time evolution of the 
yield of resonances with large decay width depends not 
only on their decay width, but also on mass difference 
between initial and final states. Similar time evolution 
occurs for the other E*, which quantitatively depends on 
their mass, decay products masses and decay width. 

For most E*s, the decay products in the channel 
A(1520) -|- TT are heavier than the decay products in oth- 
ers channels, which are thus favored by phase space. 
For most resonances in our range of temperature, the 
decay into A(1520) -I- tt remains weak. The exception 
is S(1750) which decays also to E + 77, see figure [T] 
(ms + rrir, > m-A(i52o) + w^). E(1750) begins to de- 
cay dominantly to A" (1520) at relatively low temperature 
T = 116 MeV, and continues to be produced hy T, + ij 
fusion. 

As a result, allowing all reactions to go in both direc- 
tions, the ratio A(1520)/A(1520)o has a minimum. This 
is specifically due to E(1750)) decay back to A(1520) at 
small temperatures as described above. However, when 
we satisfy Eq. ^ for dead channels the only decay occurs 
in the beginning of kinetic the dead-channel model phase. 
In that case the T^^s are smaller, and the rate of reaction 
A(1520) + TT ^ Y,* exceeds the rate for backward reac- 
tion by larger amount, compared to the scenario without 
dead channels. This amplifies the effect of A(1520) sup- 
pression. In this case, E* decay to lighter hadrons right 
after they are produced by A(1520). We can see that for 
To = 140 MeV and Tq = 160 MeV A(1520) yield is always 
decreasing in the here considered temperature range. 

For E(1385) multiplicity we find a result quite differ- 
ent from A(1520) behavior discussed here, but similar 
to what we obtained in Q by a very different method 
in a smaller basis set of states. In particular, the 
E(1385) yield is enhanced, but the maximum value of 
E(1385)/E(1385)o we find is a few percent higher, since 




140 160 
T [MeV] 

FIG. 8: The ratio E(1775)/E(1775)o as a functions of tem- 
perature T{t) for different initial hadronization temperatures 
To = 140, 160 and 180 MeV (blue/bottom, black/middle and 
red/top lines), respectively. Solid lines are for calculations 
with dead channels, dashed lines are for calculations without 
dead channels. 



we took into account the Bose enhancement of interac- 
tion rates, reaction pij) . and E* production. E(1385) 
contribution to E* production is small, compared to the 
influence of the first two effects. The time (i.e. tem- 
perature) evolution of E(1385) practically does not de- 
pend on the presence of dead channels, and the maxi- 
mum enhancement of E(1385) is even less sensitive. This 
in fact indirectly confirms that E* has a small influence 
on E(1385) multiplicity. Thus we confirm that: 

a) for To = 180 MeV E(1385) evolves with the system 
following the ambient temperature; 

b) for Tq = 160 MeV E(1385) shows some increase in 
yield; 

c) for To — 140 MeV there is a strong yield increase of 
E(1385). 

While there is little sensitivity in the yield of E(1385) 
to issue of particle momentum distribution (little dif- 
ference between the two models considered, dashed and 
solid lines), the E(1385) yield is highly sensitive to ini- 
tial hadronization condition. While for E(1385) the 
yield increases with decreased hadronization tempera- 
ture, for A(1520) the opposite is true, and in particu- 
lar the smallest final A(1520) yield corresponds to the 
smallest hadronization temperature for both models. 
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FIG. 9: The ratios E(1385)/Atot (on left) and A(1520)/Atot (on right) as a function of temperature T of final kinetic freeze-out, 
for different initial hadronization temperatures To = 140, 160 and 180 MeV (blue, black and red lines, respectively). Dashed 
lines are for calculations without dead channels, solid lines are for calculations with dead channels. The dotted purple line 
gives the expected SHM chemical equilibrium result. The dash-dot line is relative yield result result for SHM with To = T. 



C. Experimentally measurable resonance ratios 

The initial hadronization yields, which we used as a 
reference in figure [5] in order to understand the physical 
behavior, are not measurable. What is commonly used as 
a reference for the yields of single strange hyperon reso- 
nances is the overall yield of the stable A°(1115), without 
the weak decay feed from S. Aside of the initially pro- 
duced particles, the experimental yield of A°(1115) also 
includes resonances decaying during the free expansion 
after kinetic freeze-out, in particular (nearly) all decays 
of E(1385), and the experimentally inseparable yield of 
S°(1193) — > 7 -I- A° decay and the decay of any further 
hyperon resonances Y* . 

Thus we normalize our final result with the experimen- 
tally observable final Aj^^. hyperon yield: 

Atot = S°(1193) + 0.91S(1385) + A + F*. (42) 

The factor 0.91 shows that 91% of end-state i;(1385) de- 
cays to A. We also included in Atot calculations decays of 
S* ^ A + X, which makes the result slightly dependent 
on ^s/lq ratio. We use ^sHq = 1, since this ratio value 
is expected at top RHIC energy |1G |. 

As noted, A(1520) and E(1385) experimentally observ- 
able yields also include any decays which occur in the 



free-streaming post-kinetic period. Thus we have: 

I](1385)ob = E(1385) + ys*(i385): (43) 
A(1520)ob = A(1520) + yA*(i520)> (44) 

where 5^^(1385) ^(1520) hyperon multiplicities 
at kinetic freeze-out temperature, and which decay to 
11(1385) and A(1520), respectively. The multiphcities 
E(1385) and A(1520) are taken at the moment of kinetic 
freeze - out. 

In figure M we present the fractional yields 
E(1385)/Atot (left), and A(1520)/Atot (right) as a 
function of temperature of final kinetic freeze-out T. 
The results for the hadronization temperatures Tq = 140 
(blue lines), Tq = 160 (black lines) and Tq = 180 MeV 
(red lines) are shown. Solid lines are for the case with 
dead channels and dashed lines are for the case when all 
reactions are going in both directions. 

In figure[n]the green dash-dotted line is the result when 
the kinetic freeze-out temperature T coincides with the 
hadronization temperature Tq. There is no kinetic phase 
in this case, only resonances decay after hadronization. 
This result is similar to SHARE result (purple, dotted 
line). The small difference is mainly due to us taking 
into account the decays 

S(1670,1750)^ A(1520) + 7r, (45) 
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which are expected/predicted in [la |. Similarly, for 
S(1385) our results for Tq=T are different from SHARE 
results because we include the decay: 

E(1670) ^ S(1385) + 7r, (46) 

expected/predicted in [l3|- These additional resonances 
are part of current particle data set [l^l ■ 

For all initial hadronization temperatures, as the 
freeze-out temperature decreases, the suppression for 
A(1520)ob/Atot ratio is larger than for A(1520)/A(1520)o 
(at the same temperature T of final kinetic freeze- 
out). This is particularly evident for dead channels and 
hadronization temperatures Tq = 160, 180 MeV (see fig- 
urelH]). The effect is due to E(1775) suppression, as shown 
in figure[5](and similar for other S*). For Tq = 140 MeV 
the additional suppression of A(1520), described above, 
is relatively small. 

For To 140 MeV in the case without dead channels 
at final kinetic freeze-out T > 120 MeV, the final ob- 
served A(1520) suppression is even smaller, compared to 
its suppression in the kinetic phase at the same temper- 
ature (see figure [S]). The reason is that yield of E(1775) 
(and of the other S*s) is much enhanced for this range 
of temperatures see figure [8] This additional S(1775) 
decays back to A(1520). That results in a smaller sup- 
pression at these temperatures. 

The above suppression effect increases in magnitude 
for higher hadronization temperatures, since the suppres- 
sion of I](1775) and the sensitivity of A(1520)ob multi- 
plicity to Yj* decays increase with temperature. How- 
ever, when we consider dead channels (see figure [9]), the 
former effect of A(1520) suppression during evolution of 
kinetic phase increases for decreasing hadronization tem- 
peratures. Thus in the combined effect, the observable 
relative suppression of A(1520)ob/Atot, is approximately 
of the same magnitude for all hadronization temperatures 
Tq. However, the initial hadronization yield of A(1520) 
is sensitive to temperature, and decreases rapidly with 
T. Therefore only for Tq = 140 MeV, a kinetic freeze- 
out temperatures « 95 — 105 MeV, and allowing for 
dead channels, the ratio Aob(1520)/Atot reaches the ex- 
perimental domain Aob(1520)/Atot < 0.042 ± 0.01 jH, i] 
shown in figure [9] by dashed lines. 

For the same initial conditions, that is for Tq = 140 
MeV, we find the ratio E(1385)/Atot ~ 0.45 at T w 100 
MeV (and for the entire range 95 - 135 MeV, in good 
agreement with experimental data [1, 0]). In Q this 
value of S(1385)/Atot is found at T = 120 MeV, which 
was in the reference the presumed lowest possible temper- 
ature of the final kinetic freeze-out. Here we find that at 
T = 120 MeV the ratio S(1385)/Atot can be even higher 
(about 0.47), which is due to the Bose enhancement of 
in-medium S(1385) production rate (see discussion fol- 
lowing figure [5]). 



IV. CONCLUSIONS 

The resonant hadron states, considering their very 
large decay and reaction rates, can interact beyond the 
chemical and thermal freeze-out of stable particles. Thus 
the observed yield of resonances is fixed by the physi- 
cal conditions prevailing at a later breakup of the fire- 
ball matter rather than the production of non-resonantly 
interacting hadrons. Moreover, resonances, observed in 
terms of the invariant mass signature, are only visible 
when emerging from a more dilute hadron system given 
the ample potential for rescattering of decay products. 
The combination of experimental invariant mass method 
with a large resonant scattering makes the here presented 
population study of resonance kinetic freeze-out neces- 
sary. The evolution effects we find are greatly amplified 
at low hadronization temperatures where greatest degree 
of initial chemical equilibrium is present. 

Our study quantifies the expectation that in a dense 
hadron medium narrow resonances are "quenched" 
that is, effectively mixed with other states, and thus their 
observed population is reduced. Since we follow here the 
particle density, the effect we study is due to incoher- 
ent population mixing of A(1520), in particular with S*. 
This effect is possible for particle densities out of chem- 
ical non-equilibrium. However, this mixing can occur 
also at the amplitude (quantum coherent) level. As the 
result the yield suppression effect could further increase, 
in some situations further improving the agreement with 
experiment. 

Our results show that the observable ratio 
A(1520)ob/Atot can be suppressed by two effects. 
First A(1520) yield is suppressed due to excitation of 
heavy I]*s in the resonance scattering process. Moreover, 
the final A(1520)ob yield is suppressed, because I]*s, 
which decay to A(1520), are suppressed at the end of 
the kinetic phase evolution by their (asymmetric) decays 
to lower mass hadrons, especially when dead channels 
are present (see figure [5]). As a result, fewer of these 
hadrons can decay to A(1520)ob during the following 
free expansion. A contrary mechanism operates for the 
resonances such as 11(1385), A(1230). These resonances 
can be so strongly enhanced, that in essence most final 
states strange and non-strange baryons come from a 
resonance decay. 

We note that despite a scenario dependent resonance 
formation or suppression, the stable particle yields used 
in study of chemical freeze-out remain unchanged, since 
all resonances ultimately decay into the lowest "stable" 
hadron. Therefore after a description e.g. within a statis- 
tical hadronization model of the yields of stable hadrons, 
the understanding of resonance yields is a second, and 
separate task which helps to establish the consistency 
of our physical understanding of the hadron production 
process. 

We conclude noting the key result of this study, that 
we can now understand the opposite behavior of A(1520) 
(suppression in high centrality reactions) and S*^1385) 
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(enhancement, and similarly A(1230)) by considering 
their rescattering in matter. In order to explain both, the 
behavior of the A(1520)ob/Atot and E(1385)/Atot ratios, 
one has to consider T = 95—100 MeV as the favorite tem- 
perature of final kinetic freeze-out of hadron resonances, 
with To = 140 MeV being the favored chemical freeze- 
out (hadronization, QGP break-up) temperature. When 
there is little matter available to scatter, e.g. in periph- 
eral collisions, the average value of A(1520)ob/Atot ratio 
is higher, approaching the expected chemical freeze-out 



hadronization yield for Tq = 140 MeV. All these findings 
are in good agreement with available experimental data. 
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